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We study the critical and off-critical (Griffiths- McCoy) regions of the random transverse-field Ising 
spin chain by analytical and numerical methods and by phenomenological scaling considerations. 
Here we extend previous investigations to surface quantities and to the ferromagnetic phase. The 
surface magnetization of the model is shown to be related to the surviving probability of an adsorbing 
walk and several critical exponents are exactly calculated. Analyzing the structure of low energy 
excitations we present a phenomenological theory which explains both the scaling behavior at the 
critical point and the nature of Griffiths-McCoy singularities in the off-critical regions. In the 
numerical part of the work we used the free-fermion representation of the model and calculated the 
critical magnetization profiles, which are found to follow very accurately the conformal predictions 
for different boundary conditions. In the off-critical regions we demonstrated that the Griffiths- 
McCoy singularities are characterized by a single, varying exponent, the value of which is related 
through duality in the paramagnetic and ferromagnetic phases. 
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I. INTRODUCTION 



Magnetic systems with quenched disorder at very low 
or even vanishing temperature have attracted a lot of 
interest recently. In particular the quantum phase tran- 
sition occurring in quantum Ising spin glasses in a trans- 
verse field [|l| and random transverse Ising ferromagnets 
1^,^ turned out to have a number of surprising features. 
For instance the presence of quenched disorder has more 
pronounced effects on quantum phase transitions than on 
those phase transitions, which are driven by thermal fluc- 
tuations. For example in the Griffiths phase, which is at 
the disordered side of the critical point, the susceptibility 
has an essential singularity in classical systems, whereas 
in a random quantum system the corresponding singular- 
ity (Griffiths-McCoy singularity |^,||) is much stronger, 
it is in a power-law form. 

Many interesting features of random quantum systems 
can already be seen in one-dimensional models. After the 
pioneering work by McCoy and Wu |6| and later studies 
by Shankar and Murphy ||7|], Fisher [2| has recently per- 
formed an exhaustive study of the critical behavior of 
the random transverse-field Ising spin chain. He used a 
renormalization group (RG) approach, which he claims 
becomes exact at the critical point. The same type of 
method has later been used for other Id random quan- 
tum problems Q and recently some exact results are ob- 
tained through a mapping of the random XY-mode\ onto 
a Dirac equation in the continuum limit . 

In the present paper we consider the prototype of ran- 
dom quantum systems the random transverse-field Ising 
chain defined by the Hamiltonian: 



H 



I 



(1.1) 



Here the erf, trf are Pauli matrices at site I and the J; 
exchange couplings and the hi transverse-fields are ran- 
dom variables with distributio ns 7r (J) and p{h), respec- 
tively. The Hamiltonian in eq(l.l) is closely related to 
the transfer matrix of a classical two-dimensional layered 
Ising model, which was first introduced and studied by 
McCoy and Wu |]. 

In the following we briefly summarize the existing ex- 
act, conjectured and numerical r esult s on the random 
transverse- field Ising chain in eq(l.l). The quantum 
control-parameter of the model is given by 



Wnh]. 



var[ln h] 4- var[ln J] 



(1.2) 



For (5 < the system is in the ordered phase with a 
non- vanishing average magnetization, whereas the region 
(5 > corresponds to the disordered phase. There is 
a phase transition in the system at (5 = with rather 
special properties, which differs in several respects from 
the usual second-order phase transitions of pure systems. 
One of the most striking phenomena is that some physi- 
cal quantities are not self-averaging, which is due to very 
broad, logarithmic probability distributions. As a conse- 
quence the typical value (which is the value on an event 
with probability one) and the average value of such quan- 
tities is different. Thus the critical behavior of the system 
is primary determined by rare events, dominating the av- 
eraged values of various observables. 

The average surface magnetization close to the critical 
point vanishes as a power law mg ~ 5^' , where 



1 



(1.3) 



is an exact result by McCoy and Wu [||. The average 
bulk magnetization is characterized by another exponent 
/3, which is conjectured by Fisher in his RG-treatment 



(1.4) 



where r = (1 + \/5)/2 is the golden- mean. The aver- 
age spin-spin correlation function G{1) ~ [{'^i<^ij^i)]av in- 
volves the average correlation length ^, which diverges at 
the critical point as ^ The RG result by Fisher 

i is 

J^av = 2 , (1.5) 

On the other hand the typical correlations have a faster 
decay, since ^typ ~ |i5|"'''''p with vtyp = 1 0- 

In a quantum system statistics and dynamics are in- 
herently connected. Close to the critical point the re- 
laxation time tr is related to the correlation length as 
tr where z is the dynamical exponent. The ran- 

dom transverse-field Ising spin chain is very strongly 
anisotropic at the critical point, since according to the 
RG-picturc [p| and to numerical results |l0| 



Intr 



1/2 



(1.6) 



which corresponds to z = oo. On the other hand the re- 
laxation time is related to the inverse of the energy-level 
spacing at the bottom of t he s pectrum tr ~ {AE)~^. 
Then, as a consequence of eq( |l.6| ) some energy-like quan- 
tities (specific heat, bulk and surface susceptibilities, 
etc.) have an essential singularity at the critical point, 
while the correlation function of the critical energy- 
density has a stretched exponential decay, in contrast to 
the usual power law behavior. 

Leaving the critical point towards the disordered phase 
the rare events with strong correlations still play an im- 
portant role, until we reach the region S > Sg, where 
all transverse-fields are bigger than the interactions. In 
the region < S < Sq, which is called the Griffiths- 
McCoy phase the magnetization is a singular function 
of the uniform longitudinal field as rrising \Hx\^^^, 
where the dynamical exponent z varies with 6. At the 
two borders of the Griffiths-McCoy phase it behaves as 

1/26- (l + 0{S)j § as (5 ^ and z = 1 as S ^ 6^, 

respectively. 

Some of the above mentioned results have been nu- 
merically checked jl^ and in addition various probabil- 
ity distributions and scaling functions have been numeri- 
cally determined [pj]| . Our present study, which contains 
analytical and numerical investigations extends previous 
work in several respects. In the case of a limiting ran- 
dom distribution we obtain exact results on the surface 
magnetization exponent Xs — Ps/i^ and in particular the 
thermal exponent v through a simple directed walk con- 
sideration. The mapping of the problem of calculating 



various universal quantities for the random transverse 
Ising chain onto the calculation of statistical properties 
of appropriately defined random walks is one of the main 
achievements of the present paper. With its help we are 
also able to explain quantitatively many of the exotic fea- 
tures of the Griffiths-McCoy region on both sides of the 
transition. 

Moreover, we improved the accuracy of the numerical 
estimates on the bulk magnetization exponent in eq( |l.4[ ) . 
Furthermore we present new results on the magnetization 
profiles in confined critical systems as well as about the 
probability distribution of several quantities. 

Throughout the paper we use two types of random dis- 
tributions. 1) The binary distribution, in which the cou- 
plings can take two values A and 1/A with the probability 
p and q = 1—p, respectively, while the transverse-field is 
constant: 



7r( J) = p5{J - A) + qSiJ - X-^) 
p{h) = 5{h - ha) , 



(1.7) 



The critical point is given by (j) — q) In A = In ho and 
— \/vl 1 for A < 1. 2) The uniform distribution in 
which the couplings and the fields have uniform distribu- 
tions: 



7r(J) = 
p{h)^ 



1, for < J < 1 
0, otherwise 
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for < /i < ft-o 
otherwise 



(1.8) 



and the critical point is at /iq = 1 and &q = oo. 

The structure of the paper is the following. In Section 
2. we present the free-fermionic description of our model 
together with the way of calculation of several physical 
quantities in this representation. In Section 3. the sur- 
face magnetization and several critical exponents are cal- 
culated exactly using a correspondence with an adsorbing 
walk problem. Phenomenological considerations and nu- 
merical estimates about the distribution of low energy 
excitations are compared in Section 4. Numerical results 
for different critical and off-critical parameters are pre- 
sented in Sections 5. and 6., respectively. Our results are 
discussed in the final Section. 



II. FREE FERMION REPRESENTATION 



We consider the random transverse-field Ising spin 



chain in eq(l.l) on a finite chain of length L with free or 
fixed boundary conditions, i.e. with = 0. The hamil- 



tonian in eq(l.l) is mapped through a Jordan- Wigner 



transformation and a following canonical transformation 
lfl2i into a free fermion model: 



(2.1) 



2 



where 77^ and rjq are fermion creation and annihilation 
operators, respectively. The fermion energies are ob- 
tained via the solution of an eigenvalue problem, which 
necessitates the diagonalization of a 2L x 2L tridiagonal 
matrix T with non- vanishing matrix-elements T2i-i^2i — 



T2i,2i-1 — /li, « — 1, 2, . . . , L and T2i^2i+1 — 72i+l,2i — Ji, 

i = 1,2, . . . , L — 1. We denote the components of the 
eigenvectors Vq as ^,(21—1) = —4>q{i) and Vq{2i) = ijjq{i), 
i = 1,2, ... ,L, i.e. 



/ hi 
hi Ji 
Ji 



Jl-1 



Jl-1 


hL 



hL 

y 



K = 



V 



(2.2) 



One is confined to the > part of the spectrum ||Tj 



A. Magnetization 

Technically the calculation is more simple, when the 
boundary condition (b.c.) does not break the symmetry 
of the Hamiltonian, therefore we start with a chain with 
two free ends. As a consequence, in this case the ground 
state expectation value of the local magnetization opera- 
tor (Ojcrf |0)^''°° is zero for finite chains. Then the scaling 
behavior of the magnetization at the critical point is ob- 
tained from the asymptotic behavior of the (imaginary) 
time-time correlation function: 



G,{t) ^ (0|ar(r)ar(0)|0) 

= El<*l'^r|0)Pexphr(i?,~i?o 



(2.3) 



where |0) and \i) denotes the ground-state and the i-th 
excited state of H in eq(pj]), with energies Eq and Ei, 
respectively. In the thermodynamic limit in the ordered 
phase of the system the first excited state is asymptoti- 
cally degenerate with the ground-state, thus the sum in 
eq( ^.3| ) is dominated by the first term. In the large r 
limit limT-^oo G; (r) = mf, thus the local magnetization 
is given by the off-diagonal matrix-element: 



(2.4) 



In the fermion representation the magnetization operator 
is expressed as: 



erf = A1B1A2B2 . ..Ai_iBi_iAi 



(2.5) 



with 



(2.6) 



Using |1) = Vil^) the matrix-element in eq(2.4) is eval- 
uated by Wick's theorem. Since for i ^ j {0\AiAj\0) = 
{0\BiBj\0) = we obtain for the local magnetization 



m 



where 



free 



Hi Gil G12 
H2 G21 G22 



Gii- 
G21- 



Hi Gil G, 



12 



G, 



ii-i 



(2.7) 



G,k = {0\BkA,\0) 



$i(j) 

-- -J2^qik)<^qij) . 



(2., 



We note that the off-diagonal magnetization ° in 
eq(^) can be used to study the scaling behavior of the 
critical magnetization through finite size scaling. 

Next we turn to consider the system with symmetry 
breaking b.c.-s, when one of the boundary spins is fixed. 
Then one should formally put hi = ot hj^ = 0. For in- 
stance if we fix the spin at site L we put h^ = 0, implying 
that (t£ now commutes with the Hamiltonian and S'£ is 
a good quantum number. In the fermionic description 
the two-fold degeneracy of the energy levels, correspond- 
ing to iSf = -1-1 and Sf — —1, is manifested by a zero 
energy mode: £1 = in eq( p.l[ ), with the eigenvector 
Vi{i) = S2L,i, whereas for q > 1 Vq{2L) = 0. With this 
modification the ground state expectation value of the 
magnetization 



frce-t- 



(2.9) 



is formally given by the determinant in eq(2.7). The sur- 
face magnetization rus = mi = ^i(l) can be computed 
in a straightforward manner via the normalization con- 
dition J2i 'All*) — 1 leading to the exact formula [p^ : 



3 



L-1 / 



1=1 j=i ^-^j 



-1/2 



(2.10) 



where we made use of the fact that £i = because of the 
exact degeneracy of the ground state. 

Next we consider those boundary conditions, when 
both boundary spins are fixed. Then one should for- 
mally put hi — and /i^ = 0. This situation is, however, 
more complicated than the mixed b.c, since it describes 
both parallel (++) and anti-parallel (H — ) boundary con- 
ditions. To determine the magnetization profiles in these 
cases we make use the duality properties of the quantum 
Ising model. First we define the dual Pauli-operators 



^2 



4+1 



rti/2rr+i/2 , (2.11) 



in terms of which the Hamiltonian in eq(l.l) is expressed 
as: 



I 



(2.12) 



In the dual model formally the couplings and fields are 
interchanged, thus the dual Hamiltonian has zero surface 
fields, since Jq = Jl = in eq(l.l). Then it can be eas- 
ily shown that the even sector (i.e. those states which 
contain even number of fermions) of the free chain cor- 
responds to the (++) parallel boundary condition of the 
dual model, whereas the odd sector of the free chain to 
the (+-) anti-parallel boundary condition. 

To obtain the magnetization profile for fixed boundary 
spin conditions we use the expression 



'i-l-l/2 



(2.13) 



which can be obtained from eq(2.11), then express the 



product of cr^-s by fermion operators through the rela- 
tion af = AiBi and evaluate the r.h.s. of eq( p.l3 ) in 
the corresponding free-chain situation. For the parallel 
spin boundary condition we take the free-chain vacuum 
expectation value: 

m++ = (0|rf+i/2l0)++ = (0|^o5o^iSi • . .^/i3i|0)'™ , 

(2.14) 

which can be expressed through Wick's theorem as: 



1 
Gn 





Gl2 



Gil G, 
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Gil 
Gn 



(2.15) 



Here Gjk is the same as in eq(^|^), however it is calcu- 
lated with the dual couplings hi ^ Ji. 

To obtain the magnetization profile in the anti-parallel 
spin boundary condition (-1—) one s hould take the expec- 
tation value of the r.h.s. of eq( ^.13 ) in the lowest state of 



the odd sector of the free chain, which is the first excited 
state of the Hamiltonian: |1) = rji\i)). Thus: 



= (0|t;"+i/2|0)+- = {l\AoBaAiBi . . . AiBi\l 



, free 



(2.16) 

To evaluate this expression by Wick's theorem one should 
notice that the state |1) = '7]^|0) can be considered the 
vacuum state of such a system, where "qi and vii are in- 
terchanged, which formally means that ipi{k) —> —ip ijk ) 
and Ci —> — ei. Then the expectation values in eq( ^.§| ) 
are modified as 

Gjk = {l\BkAj\l) = -Y,'^,{k)<^>,{j) + ^iik)<^>iij) , 

(2.17) 

which again have to be evaluated with the dual couplings 
hi <"> Ji. Then the m'l profile is given by the deter- 



minant in eq(2.15), where the matrix-elements Gjk are 
replaced by Gjk- 



B. Susceptibility and Autocorrelations 

The local susceptibility xi at site I is defined through 
the local magnetization m/ as: 



,. 6mi 
Xi = hm — — 

Hi^i) 5 Hi 



(2.18) 



where Hi is the strength of the local longitudinal field, 
which enters in the Hamiltonian in eq(l.l) as Hiaf . xi 
can be expressed as: 



XI 



2E 



i(»krio)p 

Ei — En 



(2.19) 



which for boundary spins is simply given by: 



(2.20) 



Next we consider the dynamical correlations of the sys- 
tem as a function of the imaginary time t. First, we note 
that the correlations between surface spins can be ob- 
tained directly from eq(|2.3[) as: 



(2.21) 



Gi(r)=5]|$,(l)pexp(-re,). 



For bulk spins the matrix- element (i|crf |0) in eq(2.3) is 
more complicated to evaluate, therefore one goes back 
to the first equation of ( |2.3| ) and considers the time- 
evolution in the Heisenberg picture: 



4 



fjf (r) = exp{TH)af exp(-TiJ) K{r)<jf+n) - (Ai(T)Bi(r) • • • Ai{t)AiBi . . . Ai+^) , 

= ^i(T)Bi(r)...A,_i(r)B,_i(T)Aj(T) . (2.22) (2.23) 

can then be evaluated by Wick's theorem as a product of 
The general time and position dependent correlation two-operator expectation values, which in turn is written 
function into the compact form as a Pfafhan: 



(Ai(r)Si(T)) {Ai{t)A2{t)) {Ai{t)B2{t)) ■■■ {Ai{t)Ai{t)) (Ai(t)Ai) ••• (Ai(r)A,+„) 
(i?i(r)A2(r)) (i?i(r)i?2(T)) •.• {B,{t)Ai{t)) {B,{t)A,) ■■■ (Bi(r)A,+„) 

(A2(t)B2(t)) ••• {A2{T)AiiT)) (A2(t)Ai) ••• (A2(t)A,+„) 



±[detCy]^/^ 



(B 



l+n- 



-lA 



l+n) 



(2.24) 



where Cij is an antisymmetri c mat rix dj — —Cji, with 
the elements of the Pfaffian ( 2.24 ) above the diagonal. 
At zero temperature the elements of the Pfaffian are the 
following: 



{A,{T)Ak) = ^$g(j)$,(fc)exp(-Te,) , 
q 

{A,{T)Bk) = ^$,(j)vI/,(fc)exp(-Te,) , 

I]*9(j)*9(fc)exp(~re,) , 
q 

^*,(j)$,(fc)exp(-re,) , 



{BjiT)B,) 
{Bj{T)Ak) 



(2.25) 



whereas the equal-time contractions are given in (2.6) 
For the finite temperature contractions see c.f. p5[. 



III. SURFACE MAGNETIZATION AND THE 
MAPPING TO ADSORBING WALKS 

In this Section we analyze the surface magnetization 
of the RTIM using a mapping to an adsorbing random 
walk problem. In this way we obtain exact results for 
the critical exponents (3s and v, as well as for the surface 
magnetization scaling dimension x'^ . These observations 
will then be used in the following Section to identify the 
structure of the strongly coupled domains (SCD), which 
are responsible for the low energy excitations in the sys- 
tem. 



A. Surface magnetization and correlation length 



The surface magnetization in eq(2.10) represents per- 
haps the simplest order-parameter of the transverse-field 



Ising chain. Note that the scaling behavior of end-to-end 
correlations 



(3.1) 



is identical to that of the surface magnetization since 
Cl as well as involve only surface spin operators 
that have anomalous dimension Xs- So for instance 
Cl{5 - 0) ^ and Cl^oo{S) ~ 



i =1 2 3 4 



H — h 



H — h 



FIG. 1. Sketch of the correspondence between a bond con- 
figuration and a random walk. The thick segments in the hor- 
izontal line indicate strong (Ji — X~^) bonds, thin segments 
weak bonds {Ji = A). This bond configuration corresponds 
to a surviving walk (as indicated by the broken line staying 
completely above the horizontal line), implying a finite sur- 
face magnetization. Actually for this example ms = 1/3^^^ in 
the case A = 0, since the random walk touches the horizontal 
line three times (including the starting point). 

In the following we determine its average behavior for 
the s ymm etric binary distribution, i.e. with p=q=l/2 
in eq(L7). First we consider the system of a large, but 
finite length L at the critical p oint. In the limit X —^ 
the surface magnetization in eq(2.10), which is expressed 
as a sum of products has a simple structure. Consid- 
ering a random realization of the couplings the surface 
magnetization is zero, whenever a product of the form 
of Y\\^i "^r^' ^ ~ 1, 2, . . . , L is infinite, i.e. the number 
of A couplings exceeds the number of A^^ couplings in 
any of the [1, 1] intervals. On the other hand, if the sur- 



5 



face magnetization has a finite value, it could be of the 
form TOs = + 1)^^^, where n = 0, 1, 2, . . . measures 
the number of intervals, which are characterized by hav- 
ing the same number of A and couplings. To have 
a more transparent picture we represent the distribution 
of couplings by directed walks, which start at zero and 
make the i-th steps upwards (for a coupling Ji — X~^) 
or downwards (for a coupling Ji = X). As illustrated in 
Fig.l the surface magnetization corresponding to a ran- 
dom sequence is non-zero, only if the representing walk 
does not go below the 'time'-axis, while the correspond- 
ing surface magnetization is given as = l/(fi + 1)^^^, 
where now n just counts the number how many times the 
walk has touched the t-axis for < > 0. 

Then the ratio of walks representing a sample with fi- 
nite surface magnetization is just the survival probability 
of the walk Psurv, which is given by Ps^rviL) oc L~^/^, for 
walks (samples) of length L. In the thermodynamic limit 
this probability vanishes, thus the typical realization of 
the chain, i.e. the event with probability one has zero 
surface magnetization. This is certainly different from 
the average value, which is dominated by the rare events 
represented by surviving walks with an = 0(1). Con- 
sequently the average surface magnetization of a critical 
chain of length L is given by 



(3.2) 



On the other hand one knows from finite size scaling that 
[ms{L,S = 0)]av ^ L^"^'", where a:^ — (ig/u is the scal- 
ing dimension of the surface magnetization. Thus from 
eq(3.2) we can read the exact result: 



(3.3) 



while the prefactor in eq(3.2) is obtained as A= .6431961 
from a numerical calculation. 

In the following we argue that the exponent in eq( |3.3| ) 
is the same far all values A of the binary distribu- 
tion. First we mention that the finite-size critical surface 
magnetization [to^(L, 0)]at) is a monotonically decreas- 
ing function of A < 1: [m^^ (L, 0)]a„ < [m^^ (L, 0)]a^ < 
Ai < A2 < 1 for any value of L. Thus the corresponding 
exponents also satisfy a;^(Ai) < a;^(A2). However ac- 
cording to exact results the value of a;^ is the same, i.e. 
1/2 for the homogeneous model (A = 1) and for the 
extre me i nhomogeneous model (A — > 0), thus the relation 
in eq(3.3) should hold for any parameters of the distri- 
bution. We note that this is the first exact derivation 
of the cc^ exponent, the value of which is in agreement 
wit h pr evious ly k nown exact and conjectured results in 
eqs(L3) and (l.£). 

Next we calculate the v exponent from the 5 depen- 
dence of the surface magnetization. In the scaling limit 
L ^ 1, <C 1 the surface magnetization can be written 
as: 



Here the scaling function ms(j/), which depends on the 
ratio of L and the correlation length ^, can be expanded 
as rfis{y) = 1 + By + 0{y^), so that one obtains for the 
S correction to the surface magnetization: 



(3.5) 



with O = 1/v — xf^. This exponent can also be de- 
termined in the A — > limit of the binary distribution. 
Now, slightly outside the critical point the products of 
I terms in the sum of eq( ^.lO ) will contain a factor of 
(1 + (5)^' ~ 1 -t- 21S in leading order of S. Then the surface 
magnetization of a coupling distribution which is repre- 
sented by a surviving walk is given by: 



1 + ^(1 + 2Z,5) 



-1/2 



(l + n)-i/2_5 



{n + 1)3/2 



0(6' 



(3.6) 



where li gives the position of the i-th touching point of 
the walk with the i-axis. Next we consider a typical sur- 
viving walk, which has n = 0(1) return points, (since 
the probability of n returns decreases exponentially) and 
these points are situated at h = 0{L^^'^). Consequently 
for a typical surviving walk the correction term in eq( [3.6| ) 
is 0(L^/^), what should be multiplied by the surviv- 
ing probability of 0{L~^^'^). Since the surviving walks 
have a sharp probability distribution we are left with 
the result: K(L,(5)],,„ - [m^(L,0)]a^, = 36 + 0(6^), 
where the constant is given from numerical calculations 
as: B = 0. 270 563 ~ 17/207r. Comparing our result with 
that in eq( |3.5[ ) we get for the correlation length critical 
exponent 

v = 2. (3.7) 

This is the first exact determination of the v exponent, 
previously conjectured in Fisher's RG-study M. 



B. Relation between surface magnetization and 
adsorbing random walks 

Here we summarize and extend the mapping between 
the surface magnetization of the RTIM and the surviving 
probability of the corresponding adsorbing random walk. 
First, we consider aga in the extreme binary distribution 
of couplings in eq(1.7) with ho — 1, such that the control 
parameter in cq(1.2) is given by: 



6 = 



q-p 1 

Apq In A 



(3.8) 



(3.4) 



Then, according to the considerations of the previous sec- 
tion, the corresponding adsorbing walk has an asymmet- 
ric character: it makes steps with probabilities p and 



6 



q = 1 — p oS and towards the wall, respectively. The 
corresponding control p aram eter = q — p in eq( Al ) is 
proportional to S in eq(|3.8[). Our basic observation can 
be summarized as: 



ms{5,L) 



J 



(3.9) 



At the critical point from eq(A7) Psurv(0,i) ~ L^'^ 
and "f — 1/2 is just a;f„, the surface magnetization scal- 
ing dimension of the RTIM. In the paramagnetic phase of 
the RTIM S > and the corresponding walk has a drift 
towards the adsorbing wall. The surviving probability in 



^'surv(<5^ > 0,L) exp{-L/^^) 

e 



r-2 

2 ' "w 5 



(3.10) 



is characterized by a correlation length, which diverges 
as ^ with v = 2. We note that the expression for 
iy, in ( |3.10D agrees with the RTIM result by Fisher 1). 
Finally, in the ferromagnetic phase 5 < Q the correspond- 
ing walk drifts off the wall and the surviving probability 
has a finite limit as L ^ oo: 



PsurviS < 0,L-> oo) 



p-q 

P 



(3.11) 



This expression then corresponds to a finite average sur- 
face magnetization of the RTIM, which linearly vanishes 
at the critical point. Thus the surface magnetization ex- 
ponent of the RTIM is /3s = 1, in agreement with the 
exact results by McCoy and Wu 

These results obtained for the extreme binary distribu- 
tion can be generalized for other random distributions, 
too. It is enough to notice, that the surface magneti- 
zation in a sample with non-surviving walk character 
is exponentially vanishing with the size of the system. 
Therefore the basic relation in eq(3.9) remains valid. 
Then at the critical point the Psuiv('5 = 0, i) ~ L^^/^ 
is a consequence of the Gaussian nature of the random 
walk. Similarly, the relations ~ > and 

^surv(^tu < 0, L) ~ —Sw, (5u, < follow from the scahng 
behavior of the random walks. 



IV. DISTRIBUTION OF THE LOW-ENERGY 
EXCITATIONS 

In the previous Section we saw that the surface order 
is connected to such a coupling distribution, which can 
be represented by a surviving walk. Since (local) order 
and small (vanishing) excitation energies are always con- 
nected, we can thus identify the local distribution of cou- 
plings, which result in a strongly coupled domain (SCD). 
Note that a SCD is not simply a domain of strong bonds, 
but it generally has a much larger spatial extent. 

To estimate the excitation energy e of an SCD we make 
use the exact result for the lowest gap €i{l) of an Ising 



quantum chain of I spins with free b.c. p^ : Since we are 
interested in a bond- and field configuration that gives 
rise to an exponentially small gap ei we can neglect the 
r.h.s. of the eigenvalue equation 



(4.1) 



c.f. eq(2.2) and derive approximate expressions for the 
eigenfunctions (f>i and ^'i. With these one arrives at 



ei(0 



(4.2) 



Here rUg and mj denote the finite-size surface magneti- 
zations at both ends of the chain, as defined in eq( 2.1C ) 
(for toJ simply replace hj/Jj by Hl-j / Jl-j in this eq.). 
If the sample has a low-energy excitation, then both ends 
surface magnetizations are of 0(1), consequently the cou- 
pling distribution follows a surviving wal k picture. The 
corresponding gap estimated from eq( |4.2| ) is given by: 



ei 



/-I 

n 



J,. 



exp 



{-/tr-ln(J//i)} , (4.3) 



where Itr measures the size of transverse fluctuations of 
a surviving walk of length I and \n{J/h) is an average 
coupling. In the following we assume that the excitation 
energy of surface SCD-s are of the same order of mag- 
nitude as those localized in the bulk of the system and 
have the same type of coupling distributions. Thus we 
identify the SCD-s, both at the surface and in the vol- 
ume of the system, as a realization of couplings and fields 
with surviving walk character and having an excitation 
energy given eq(4.3). With this prerequisite we are now 
ready to apply our theory for the critical and off-critical 
regions of the RTIM of L sites. 

At the bulk critical point the characteristic length / of 
surviving regions is of the order of the size of the system 
L, thus the SCD extends over the volume of the system. 
The transve rse fl uctuations of the couplings in the SCD 
are from eq(Al3) as Itr ~ L^^"^, thus we obtain for the 
scaling relation of the energy gap at the critical point: 



e{6 = 0, L) - exp(-const • L^/^) 



(4.4) 



in accordance with the existing numerical results pC| ] . At 
this point it is useful to point out the origin of the expo- 
nent 1/2 accompanying the length scale L in eq(4^): it is 
the fact that the sequence of hi / Ji is random and uncor- 
related, for a general sequence one would have Itr ~ L'^ 
with Lj being the wandering exponent (the scaling di- 
mension of the transverse fluctuations) of the particular 
sequence under considerations. For instance one could 
also consider relevant aperiodic sequence (generated in a 
deterministic fashion), which have either the same or dif- 
ferent wandering exponents, leading either to the same 
or a different scaling behavior of the energy scale at the 
critical point as the random chain studied here |18[. 
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FIG. 2. The integrated gap probability distribution 
r2_L(lnei) in the disordered phase {h > 1) for different val- 
ues of h. The dynamical exponent zijn) is extracted from the 
expected asymptotic form lnr2L(lnei) = l/2:(/i) Inei) -I- const 
which is a straight line when using a logarithmic scale on the 
y-axes. Thus \lz(K) ^ 0.82, 0.62 and 0.40 for h = 4.0, 2.0 
and 1.5, respectively. The data for the uniform distribution 
averaged over 50000 samples. Note that for large h, i.e. far 
away from the critical point, there is essentially no system size 
dependence, whereas closer to the critical point the asymp- 
totic slope can only is reached only for large enough system 
sizes. 



In the paramagnetic phase the probability to find a 
SCD of size I, which is localized at a given point is pro- 
portional to exp(— Z/^), c.f. cq( |3.10 ). Since the SCD can 
be located at any point of the chain, the actual prob- 
ability is proportional to the length of the chain, thus 
Pl{1) ~ iexp(— //^). The characteristic size of SCD-s 
obtained from the condition Pl{1) = 0{1) is given by: 



I - ^Ini , S>0 



(4.5) 



which grows very slowly with the linear size of the sys- 
tem. The characteristic transv erse f luctuations of such a 
walk is given - according to eq( A16D - as Itr ~ (g — p )la, 
with a = 0(1). Putting this expression into eq( |4.3| ) we 
obtain the scaling relation 



(4.6) 



where the dynamical exponent z{S) = 2a/S is a contin- 
uous function of the control parameter S. Our estimate 
qualitatively agrees with Fisher's result ||], that close to 
the critical point z{S) = 1/2(5. 

The dynamical exponent z{d) is conveniently measured 
from the scaling behavior of the probability distribution 
Pl{1) ^ ^L(ln(e)) ~ L. For a given large L the scaling 
combination from eq([4.6|) is Le^/^, thus 



Pie) 



(4.7) 



The distribution function of the gap in eq(4.7) has al- 
ready been studied in jl^ for periodic boundary condi- 
tions. Here we considered free chains and investigated 
the accumulated probability distribution 



r!L(lne) 



dyPUy) 



(4.8) 



As seen on Fig. ^ the accumulated probability distribu- 
tion for low energies is approximately a straight line on a 
log-log plot and from the slope one can estimate l/z{S) 
quite accurately. 
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FIG. 3. The same as in the Fig. 2 in the ordered phase 
{h < 1). The data do not scale with Inei/L, there are strong 
logarithmic corrections. A scaling with Inei /[L/ln{L)] is also 
poor (as can be seen in the figure), most probably higher 
powers of ln(Z/) are involved. 
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In the ferromagnetic phase of the RTIM the size of the 
SCD is of the order of the sample, I ^ L and also the 
transverse fluctuations of the couplings are In ~ L. Con- 
sequently the energy of the first excitations scale expo- 
nentially with the size of the system: e ~ exp(— constL). 
Here, however, one should take into account that - due 
to the duality relation in eq( 2.12| ) - in a strongly coupled 
environment there are always weakly coupled domains 
(WCD), which are the counterparts of the SCD in the 
paramagnetic phase. The characteristic size of a WCD 
is O(lnL), and their presence will reduce the size of the 
SCD-s, such that one expects logarithmic corrections to 
the size of transverse fluctuations ^tr- Indeed the numer- 
ical results on the accumulated gap distribution function 
in Fig. m can be interpreted with the presence of such 
corrections. 
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3 for the second 
One observes that 



FIG. 4. The same as in Fig, 
lowest excitation, i.e. r2_L(hie2). 
asymptotically lnf2_L(lne2) = 1/ z{h)\ne2 4- const, with 
z(h = 0.5) = 0.62 — z(h — 2.0) as one would expect from 
duality, by which z{h) = z{l/h). 

In the ferromagnetic phase many physical quantities 
(connected autocorrelation function, susceptibility, etc.) 
are connected with the distribution of the second gap. 
Unfortunately, we can not make an estimate for 62 on 
the base of our present approach. However, our model is 
self-dual, the distributions of t he co uplings and the fields 
transform to each other in eq(2.12) for S —5. There- 
fore, we assume that the scaling behavior of ei in the 
paramagnetic phase and that of €2 in the ferromagnetic 
phase are also related through duality, thus 



62(6 < 0,L) ~ L 



-z{-5) 



(4.9) 



Indeed, as seen on Fig. ^ the scaling relation in eq(|4.9|) 

is satisfied, however with strong, logarithmic corrections. 



V. CRITICAL PROPERTIES 

A. Surface magnetization - canonical vs. 
micro-canonical ensemble 

The surface magnetization of the RTIM has already 
been studied in Section 3. Here we revisit this problem 
in order to answer to the question, whether the values 
of the average quantities and the corresponding critical 
exponents depend or not on the ensemble used in the 
calculations. Our present study is motivated by a recent 
work in which finite-size scaling methods and their 
predictions for critical exponents pO| ] have been scruti- 
nized for random systems. 

In our approach in Section 3 the bond- and field con- 
figurations were taken completely random according to 
the corresponding distribution. We call this the canoni- 
cal ensemble, since only the ensemble average of In Ji and 
In hi is held fixed. One can also confine oneself on a sub- 
set of this ensemble, in which the the product of all bonds 
in the chain is exactly equal to the product of all fields 
in the chain, i.e. ^f^i In Ji = X^tTi^ (note that we 
study the surface magnetization in a chain with one fixed 
boundary condition such that there are exactly as many 
bonds as fields). This wc call the micro- canonical ensem- 
ble. The motivation for the introduction of this ensemble 
can be found in essentially it is a more restrictive 
way of fulfilling the criterium [In J]av = [ln/i]av for being 
at the critical point. 

The c ritic al exponents of the canonical ensemble 
eqs(^r^),(1.5), which agree with other exact and RG re- 
sults, are then the canonical ones. 

In the following we c a lcula te the critical exponents cor- 
responding to cqs(3.3),( |l.5| ) also in the micro-canonical 
ensemble. We use again the symmetric binary distribu- 
tion in eq(1.7), such that the samples have the same num- 
ber of A and couplings. In the extreme limit A ^ 0, 
as in the canonical case, the critical point surface magne- 
tization can be determined exactly, through studying the 
surviving probability of the corresponding adsorbing ran- 
dom walk. To determine the micro-canonical surviving 
probability first we note that from the canonical walks 
only a fraction of 0{L~^^^) is micro-canonical. Second, 
the micro-canonical surviving walks have their end at 
the starting point. Such returning surviving walks are 
of a fraction of 0{L~^^'^) among the canonical walks. 
Thus the surviving probability of micro-canonical walks 
is Psurv(m.c.) ^ L^^, therefore the micro-canonical sur- 
face magnetization scaling dimension is given by: 



,(m.c.) = 1 . 



(5.1) 



The correlation length critical exponent v i s again ob- 
tained by analysing the expression in eq( |3.6| ). The typ- 
ical number of return points of the surviving walks is 
again n = 0(1), but now U = 0{L), since the end-point 
of the walk is a retu rn point. Consequently, the correc- 
tion term in eq(|3.6[) for surviving walks is 0{L), what 
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should be multiplied by the surviving probability to ob- 
tain the average of 0(1), from which the micro-canonical 
correlation length exponent 



z/(m.c. 



1 



(5.2) 



follows. Compari ng t h e ca nonical a nd m i cro- canonical 
exponents in eqs(|3.3[),(|l.5D and eqs(5J_),(5^), respec- 
tively, we can conclude that they are different. It is im- 
portant, however, to note that the surface magnetization 
exponent, defined via 



(5.3) 



is the same for both ensemles: (3s = x'mi' = 1. In Fig. |^ 
we show the scaling plots for the surface magnetization 
in the two ensembles obtained numerically by evaluat- 
ing eq( 2.1C ) for the binary distribution. Note that we 
expect similar results to hold for end-to-end correlations 
[(cr^CT£)]av, with the exponent Xg replaced by 2xs (c.f. 
sec. 3). 
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FIG. 5. Scaling plot of the surface magnetization in the 
canonical ensemble (top) and the micro-canonical (bottom) 
ensemble. Both data are for the uniform distribution aver- 
aged over 100000 samples. 

The reason for the difference of the critical exponents 
measured in the two ensembles is the fact that several 



physical quantities, among those the surface magnetiza- 
tion, is not self-averaging at the critical point. To illus- 
trate this property in Fig. || we have plotted the prob- 
ability distribution of the surface magnetization in the 
two ensemles. Both scale, as expected, as 



on ^ 1 ~ / In 
P^lnr..)^— 



(5.4) 



the asymptotic form of the scaling function p for the 
canonical ensemple was determined analytically for 
particular distributions of the fields and/or couplings it 
can even be calculated exactly . 
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FIG. 6. Scaling plot of the probability distribution 
P(lnms) of the surface magnetization in the canonical ensem- 
ble (top) and the micro-canonical (bottom) ensemble. Both 
data are for the uniform distribution averaged over 500000 
samples. 

The average is determined by the rare events having 
a magnetization of order 0(1), i.e. by the asymptotic 
behavior of p(y) for y — > 0. From fig. ^ (top) we con- 
culde that for the canonical ensemble p(?/) approaches a 
constant for y — s- 0_, whereas for the micro-canonical en- 
semble in fig. ^ (bottom) the scaling functions shows a 
power law dependence 



for 



(5.5) 
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with a positve exponent a. For the average then follows 



m, 



,]av = -^y" dm p(ln ms/\/L) (5.6) 

^-^J dm{\ lnm,|/%/I)" oc i-(i+'^)/2 

As we said above a — for the canonical ensemble, result- 
ing in [rn^lav ~ and a = 1 fits the data reasonably 
well in case of the micro-canonical ensemble, resulting in 

Based on our observation on the surface magnetization 
we assume that for other non self-averaging quantities the 
corresponding critical exponents are generally ensemble 
dependent. In the rest of the paper we restrict ourselves 
to the canonical ensemble. 



B. Profiles of observables 



where the scaling function Gab{l/L) depends on the uni- 
versality class of the model and on the type of the bound- 
ary condition. With symmetric boundary conditions the 
scaling function is constant Gaa — A. For conformally 
invariant, non- symmetric boundary conditions the scal- 
ing function has been predicted for several models. For 
the Ising model the magnetization profiles with free-fixed 
(f-|-) and (-1--) boundary conditions are predicted as: 



G 



f+ 



A 



2L 



(5.10) 



and 



Acqs ■ 



(5.11) 



A real system is always geometrically constrained and 
due to modified surface couplings its properties in the 
surface region are generally different from those in the 
bulk. Close to the critical point this surface region, which 
has a characteristic size of the correlation length, intrudes 
far into the system. At the very critical point the appro- 
priate way to describe the position dependent physical 
quantities is to use density profiles rather then bulk and 
surface observables. For a number of universality classes 
much is known about the spatially inhomogeneous be- 
havior, in particular in two-dimensions, where conformal 
invariance provides a powerful tool to study various ge- 
ometries p2| . 

In a critical system confined between two parallel 
plates, which are at a large, but finite distance L apart, 
the local densities $(r) such as the order parameter 
(magnetization) or energy density vary with the distance 
I from one of the plates as a smooth function oil/ L. Ac- 
cording to the scaling theory by Fisher and dc Gennes 



(5.7) 



where is the bulk scaling dimension of the operator 
$, while ab denotes the boundary conditions at the two 
pla tes . In a d-dimensional system the scaling function in 
eq(^?^) has the asymptotic behavior: 



Fab{l/L) = A 



l + B„ 



I 



< 1 . (5.8) 



Here the amplitude of the first correction term is uni- 
versal, the corresponding exponent is just = d the 
surface scaling dimension of $. 

In two-dimensions conformal invariance gives further 
predictions on the profile: 



L . I 

— smTT— 

TT L 



Gab{l/L) 



(5.9) 
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FIG. 7. Scaling plots of the magnetization profiles for 
non-symmetric boundary conditions {I' = I — 0.5). Top : 
plus- minu s (-1-—) b.c, the broken line is a fit to the form (5.!:) 
and (5.11) with A as a fit parameter. Bo ttom : free - fixed (f+) 
b.c, the broken line is a fit to the form (5.9) and ( |5.1C| ) with 
A as a fit parameter. 
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In two-dimensions conformal invariance can also be 
used to predict the critical off-diagonal matrix-element 
profiles (<i>|$(^)|0), where ($1 denotes the lowest ex- 
cited s tate leading to a non-vanishing matrix-element 
(see eq([2.4|)). These ofF-diagonal profiles give information 
about the surface and bulk critical behavior via finite-size 
scaling while avoiding the contribution of regular terms. 
With symmetric boundary conditions one obtains for the 
profile pi: 



($|<i>(/)|0) « (-) 



sin 4 



(5.12) 



which involves both the bulk and surface scaling dimen- 
sions. 
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FIG. 8. Top: scaling plot of the magnetization profile for 
symmetric (here; fi xed) boundary conditions. The broken line 
is a fit to the form (5.9) and Gaa = A with A as a fit param- 
eter. Bottom: the profile of the off-diagonal matr ix ele ment 
with free b.c. The broken line is a fit to the form (5.12). 



The numerically calculated average diagonal and off- 
diagonal magnetization profiles (see sec. 2) for the RTIM 
are pre sented in Fig. |^ and ^for the uniform distribution 
(in p5[ some data for the binary distribution have been 
presented). Here we do not and a;^ as fit param- 

eters, but fix them to the theoretically predicted values 



in eqs(1.4), (3.3). The only fit parameter is the non- 
universal prefactor A, which is found remarkably con- 
stant for the different boundary conditions. As one can 
see on Fig. Q the data for different length L collapse to 
scaling curves, which are very well described by the scal- 
ing functions predicted by conformal invariance. Thus 
we can conclude that not only th e sc aling prediction by 
Fisher and de Gennes |2^ in eq(5/7) is very well satis- 
fied for the RTIM, but the corrections to the appropriate 
conformal results are also very small, practically negligi- 
ble. This is an unexpected result, since the RTIM is not 
conformally invariant, due to anisotropic scaling at the 
critical point eq(1.6). 

To close this section we present numerical results for 
the bulk magnetization scaling dimension Xm and com- 
pare i t w ith Fisher's perhaps most striking prediction 
in eq(L4). Here we have made effort to increase the 
numerical accuracy, th erefo re we worked with the bi- 
nary distribution in eq(l/7) on chains with both ends 
fixed with length L < 24 and performed the exact av- 
erage of the local magnetization on the central spin. 
From the finite-lattice magnetizations, which scales as 
[m(L, 0)]av ~ Ij-^rn. -^g have determined x^ by two- 
point fit, comparing systems with sizes L and L — 2. 
From the finite-size exponents, presented in Table 1 for 
A = 2,3 and 4 one conclude, that they are in agreement 
with Fisher's result: = (3 - V5)/4 = .191. Unfortu- 
natelly the numerical data in Table 1 show log-periodic 
oscillations, which is a consequence of the energy scale in- 
troduced by the binary distribution. Therefore one can 
not use the accurate sequence-extrapolation methods to 
analyse the limiting behavior of the series. Then from a 
simple linear fit one can obtain the estimate: 

a;„, = 0.190 ± 0.003 , (5.13) 

improving the accuracy of previous MC estimates |10[. 
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.199992 


22 


.181044 
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TABLE I: Numerical estimates for the bulk magnetiza- 
tion exponent Xm{L) for the binary distribution for var- 
ious values of A. 



C. Dynamical correlations 

Th e gen eral time and position dependent correlations 
in eq(2.23) have a complicated structure at the critical 
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point. Therefore we consider dynamical correlations on 
the same spin, which has a simpler asymptotic behavior. 
First we consider the bulk autocorrelation function 



G(r) = WLi^{T)ali^) 



(5.14) 



and recapitulate the scaling argument in Ref ||26|| . 

The autocorrelation function, like to the (local) mag- 
netization, is not self-averaging at the critical point: its 
average value is determined by the rare events, which 
occour with a probability Pr and Pr vanishes in the ther- 
modynamic limit. In the random quantum systems the 
disorder is strictly correlated along the time axis, conse- 
quently in the rare events with a local order, i.e. with 
a finite magnetization also the autocorrelations are non- 
vanishing. Under a scaling transformation, when lengths 
are rescaled as I' = l/b, with 6 > 1 the probability of the 
rare events transforms as P^ = 6"^™, like to the local 
magnetization. As we said above the same is true for the 
autocorrelation function: 

G(lnr) = r^"G(lnT/&i/2) S^O, (5.15) 

where we have made use of the relation between rele- 



vant time tr and length ^ at the critical point in eq(1.6) 
Taking now the length scale as 6 = (Inr)^ we obtain: 



G(r)^ (In r) 



-2x„ 



6 = 



(5.16) 



For surface spins in eqs( |3.15 , 5.16| ) the surface magneti- 
zation scaling dimension appears. 
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FIG. 9. Bulk spin-spin autocorrelation function G(r) 
eq( 5.14 ) for various system sizes (and the uniform distribu- 
tion). The straight line is the prediction according to eq( 5.1(3) 



straig ht li ne in the infinite system size limit according 
to eq( ^.l(j ) . As can be seen in Fig. ^ the data agree well 
with this prediction. For the surface autocorrel ation s 
Gi(t) = [(erf (T)cr^)]av, evaluated according to eq( 2.21 ), 
which is much less involved than the computation of a 
pfaffian, a similar plot with the bulk magnetization expo- 
nent Xm replaced by the surface magnetization exponent 
gives also an excellent agreement with the prediction 
Gi(lnT) - (lnT)-2x^. 

To complete the results on critical dynamics we men- 
tion the scaling behavior of the energy-energy autocorre- 
lation function [{o'f {T)af)]av As shown in Ref this 
quantity is self-averaging at the critical point and can 
be characterised by a power law asymptotic decay with 
novel critical exponents, which are different in the bulk 
and at the surface of the system. 



VI. OFF-CRITICAL PROPERTIES 



A surprising property of random quantum systems is 
the existence of Griffiths-McCoy singularities in the para- 
magnetic side of the critical point. In the corresponding 
Griffiths-McCoy region the autocorrelation function de- 



cays as a power G(r) t' 



-1/^(5) 



where the dynamical 



exponent z{S) character ises also the distribution of low 
energy excitations in eq(4.7). As a consequence, the free 
energy is non-analytic function of the magnetic field and 
the susceptibility diverges in the whole region. 

According to the phenomenological theory [|l^ in the 
Grifhths-McCoy region the singularities of all physical 
quantities are entirely characterised by the dynamical ex- 
ponent z{6). Numerical calculations [ p^|jTl| ] give support 
to this assumption, although there are discrepancies be- 
tween the values of z{S) obtained from different quanti- 
ties. 

Here we extend previous investigations in several re- 
spects. First, we consider also the surface properties, 
such as the surface autocorrelation function and the sur- 
face susceptibility. Second, we investigate also the ferro- 
magnetic side of the criticalpoint. In the neighborhood 
of the critical point Fisher [|| has already obtained some 
RG results in the ferromagnetic phase. Here we are going 
to check these results numerically and to extend them for 
finite (5 < 0. 



A. Phenomenological scaling considerations 



In Fig. y we present the numerical results for the 
critical bulk autocorrel ation function obtained via eval- 
uating the Pfaffian eq( ^.24| ). Note that we have cho- 
sen L to be odd, so that L/2 denotes the central spin, 
representing the bulk behavior in a system with free 
b.c. A plot of G(r)~^/^'^", with Xm as in (5.13), ver- 
sus Inr (or r on a logarithmic scale) should yield a 



As already shown in Section 4 the dynamical exponent 
z{S) is conveniently measured from the probability distri- 
bution of the energy gap (in the ferromagnetic phase one 
considers the second gap in a finite system, which does 
not vanish exponentially.) For large systems the gap- 



distribution is given by eq(4.7) and with this the average 
autocorrelation function is given by: 
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G(7 



P(e) exp(— Te)(ie ~ t 



(6.1) 



In a finite system of size L for long enough time, such 
that T 3> L^^^\ the decay in eq(6.1) will change to a 
G{t) ^ 1/t form, which is characteristic for isolated 
spins. It means that in the above limit the system can 
be considered as an effective single spin. 

In the following we present a simple scaling theory 
which explains the form of the asymptotic decay in 
eq(|6?l|). Here in the Grif fiths- McCoy phase we modify 
the scaling relation in eg ( ^.15 ) by two respects. First, 
the scaling combination is changed to r/h^ , since the dy- 
namical exponent z(5) is finite in the off-critical region. 
Second, the rare events^ which are responsible for the 
GrifRths-McCoy singularities are now samples with very 
low energy gaps and their number is practically indepen- 
dent of the size of the system. Consequently the rescaling 
prefactor is b^^ and the scaling relation is given by: 



G{T,l/L)^b-^G{T/h\h/L) 



(6.2) 



where the inverse size of the system 1 /L is also included 
as a scaling field. Now taking b — t^/^ we obtain: 



G(r,l/L) = r-i/^G(ri/7L) 5^0, 



(6.3) 



thus in the thermodynamic limit we recover the power 
law decay in eq(|6.lj). The scaling function G{y) in eq(|6.3|) 
should behave as G{y) ~ for large y, in this wa y on e 
recovers the limiting 1/t decay, as argued below eq(3.1). 
Then the finite size scaling behavior of the autocorrela- 
tion function is of the form of L^~^ , and after integrating 
G(r, 1/ L) by r the same scaling behavior will appear in 
the local susceptibility: 



5^0 



(6.4) 



B. Numerical calculation of the dynamical exponent 



The phenomenological description of the Griffiths 
phase suggests that all Griffiths-McCoy singularities 
emerging in temperature, energy, time or frequency de- 
pendent quantities should be paramtrizable by a single 
dynamical expontent z{5). In this subsection we present 
the results on our numerical estimates for z{5) resulting 
from the calculation of the following quantities: 



• distribution of low energy excitations, 

• autocorrelation function on bulk and surface spins, 

• distribution of surface susceptibilities. 
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FIG. 10. Integrated probability distribution of the zero fre- 
quency surface susceptibility for different system sizes in the 
disordered phase at /i = 2.0. Note that z{h), as determined 
from the slope of the straight line, turns out to be within the 
error margin of z{h) determined via the gap distribution (see 
Fig. I). 
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FIG. 11. The bulk autocorrela- 

tion function [(cr2/2(''')''"L/2('^))]'»^ imaginary time in the 
disorder ed p hase {h > 1), calculated at a central spin i = L/2 
with eq(2.2c) via the PfafHan method presented in section II. 
The straight lines are fits to the expected power law decay 
.^-i/^Ci) 
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The distribution functions for the energy gaps have al- 
ready been presented in Section 4. T he same quantity for 
the surface susceptibihty in eg ( ^.20 ) has a similar form 
as the inverse gap, as seen in Fig. [l^. The only dif- 
ference that for the su scept ibility the matrix-element in 
the denominator of eq( ^.20 ) select one special position of 
the SCD. As a consequence the corresponding probability 
distribution has no L-dependence, as already discussed in 
pot and can also be seen by comparing eq(4.6) with (3.4). 
The z{S) exponents calculated from the surface suscep- 
tibility distribution agree well with those obtained from 
the gap distribution. 
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FIG. 12. The surface autocorrelation function 
[(erf (r)af (0))]av in imagi nary time in the disordered phase 
{h > 1), calculated via eq(2.21). The straight lines are fits to 
the expected power law decay t~ 



-l/z(h) 



The average autocorrelation function is measured at 
two sites of the chain: on the central spin, giving an es- 
timate for the bulk correlation function and on the sur- 
face spin. The average bulk autocorrelation functions 
are drawn on a log-log plot in Fig. |ll| for several val- 
ues of 5 > 0. One can easily notice an extended region 
of the curves, which are well approximated by straigth 
lines, the slope of which is connected to the dynamical 
exponent through eq(^^). Similar behavior can be seen 
on Fig. O where the average surface autocorrelation 



functions are drawn. Our investigation on the dynamical 
exponent is completed by studying the connected surface 
autocorrelation function in the ferro mag netic phase. As 
seen on Fig. |l3| the scaling form in eq( 11 ) is well satisfied 
for this function, too. 
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FIG. 13. The connected part of the surface autocorrela- 
tion function [({af(r)af(0))>]av = [{af{r)af{0))-mlU in 
im agina ry time in the ordered phase {h < 1), calculated via 
eq( 2.21 ), but now substracting |$i(l)p. The straight lines 
are fits to the expected power law decay t-~i/^(''). 
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FIG. 14. Summary of all estimates of the dynamical expo- 
nent z{h) as a function of the distance from the critical point 
5 — \ \n h\/2. The open circles are the estimates from the gap 
distribution, as exemplified in Fig. H, open tiangles: surface 
susceptibility distribution (c.f. Fig. Im, full circles: surface 
autocorrelation function (c.f. Fig. [12), crosses: surface auto- 
correlation function in the ordered phase (c.f. Fig. H ), open 
squares: bulk autocorrelation function (c.f Fig. [L1[) . Note 
that wheras all former estimates agree within the error mar- 
gin (which is roughly the size of the symbols) the latter esti- 
mate, namely the one obtained via the bulk autocorrelation 
function, differs significantly from all others. 



The behavior of the dynamical exponents calculated 
by different methods are summarized in Fig. |l^. First 
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we note that the numerical estimates are very close to 
each other. The only exception is the data obtained 
from the bulk autocorrelations. To explain the possi- 
ble origin of this discrepancy we turn in the Discussion. 
The z{6) values well satisfy the two theoretical limits: 



limi 



<5|- 



,z{S) = 1 and lim|5|^o = 1/2(5 



Further- 



more the dynamical exponents show the duality relation: 
z{6) = z{-6). 



VII. DISCUSSION 



and in the volume o f the sy stem. Then, from the argu- 
ment leading to eqs(15 4~6| ) one can obtain logarithmic 
corrections between the dynamical exponents. This fact 
can then explain the differences in the finite-size data. 
We have here, at the first time, numerically studied the 
GrifBths-McCoy singularities in the ferromagnetic phase, 
too. In this region the second gap of the Hamiltonian and 
the connected autocorrelation function scale with the dy- 
namical exponent, which, according to numerical results, 
satisfies the duality relation. 



In this paper the critical and off-critical properties of 
the random transverse-field Ising spin chain are stud- 
ied by analytical and numerical methods and by phe- 
nomcnological scaling theory. The previously known ex- 
act §, RG § and numerical results [|o|j25|j2|,|lll about 
the model have been extended and completed here in 
several directions. The scaling behavior of the surface 
magnetization is obtained through a mapping to an ad- 
sorbing random walk and the critical exponents /S^, v and 
xf^ are calculated exactly. We have also shown that the 
values of these exponents depend on the ensemble over 
which the averaging is performed. 

Using the correspondence between the surface mag- 
netization and the adsorbing walks we have identified 
strongly coupled domains in the system, where the cou- 
plings have a surviving walk character, and estimated 
the distribution of low energy excitations both in the 
critical and off-critical regions. This provides a com- 
prehensive explanation of the microscopic origin of the 
Griffiths-McCoy singularities. It turns out that most of 
the astonishing features of the critical as well as the off- 
critical (Grifliths-McCoy) properties can be simply ex- 
plained via random walk analogies. However, one pre- 
diction by Fisher |^ , namely the exact value of the bulk 
magnetization exponent P and its surprising relation to 
the golden mean, still lacks a simple explanation in terms 
of universal properties of random walks. 

In the numerical part of our work we have treated rela- 
tively large {L < 128) finite systems. At the critical point 
we have calculated the magnetization profiles for differ- 
ent boundary conditions, which are found to follow accu- 
rately the conformal predictions, although the system is 
not conformally invariant. We have also increased the nu- 
merical accuracy in the calculation of the bulk magnetiza- 
tion scaling dimension. In the off-critical regions we have 
determined the dynamical exponent z{d) from different 
physical quantities. The obtained results give support to 
the scaling prediction that the Griffiths-McCoy singular- 
ities are characterised by the single parameter z{d). Here 
we note that the numerical data show systematic differ- 
encies, when z(S) is calculated from bulk or from sur- 
face quantities. Similar observation has been made in ref 
PH , too. The possible origin of the discrepancies is, that 
the SCD-s, which are responsible for the GrifRths-McCoy 
singularities, have different environments at the surface 
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APPENDIX: ADSORBING RANDOM WALKS 

Here we summarize the basic properties of one- 
dimensional random walks in the presence of an ad- 
sorbing wall. For simplicity, first we consider a walker, 
which makes steps of unit lengths with probabilities p 
and q = 1 — p to the positive and to the negative direc- 
tions, respectively. Starting at a distance s > from an 
adsorbing wall we are interested in the surviving proba- 
bility Psuivi^w, L) after L steps. Here 



(Al) 



measures the average drift of the walk in one step: for 
(5^ < ((5 > 0) the walk has a drift towards (off) the 
wall. 

The probability Wl{1), that the walker after L steps is 
at a position I < L, can be easily obtained by the mirror 
method l|27|l: 



Wl{1) 

L + l L 

p 2 q 



(A2) 



LI 



(^)!(^)! 



' L+l 
V 2 



s)!( 



L-l 
2 



In the following we take s = 1 and in the limit L ^ 1, 
Z 3> 1 we use the central limit theorem to write in the 
continuum approximation / x and Wl{1) Pl{x) as: 
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r . la; 1 
Pl(x) = --—=== exp 
P L V27rLo-2 



with 



(a; — x)"^ 



(A3) 



X = [p — q)L = -~6wL, = ipq . 

The surviving probabihty is then given by: 



(A4) 



J 



f°° I 20-2 

,((5u,,L) = J PL{x)dx ^ J -^^^^e-K-p{' 



(A5) 



where 



y = 



V2Ict2 



= -<5„ 



L 

2^ ' 



(A6) 



and (jjz = 2/y^ JJ^ exp(— i2)c?t is the error function 
In the following we evaluate Psm v {5w , L) in eq([A5|) in 

the different limits. 

In the symmetric case 5^, = 0, p = 1/2: 

PsnASu, = 0, L) = -1= - _ 

For (5^ > 0, when p < q and the walk has a drift towards 
the wall the surviving probability has an exponential de- 
cay as y ^ — oo: 



PsuiviSw > 0,L) = 



2a2 _ 1 

L-i/2exp(-L/e^)^,,/L 



(A8) 



with a correlation length: 



(A9) 



Finally, for 6^ < 0, when p < q and the walk is drifted 
from the wall the surviving probability has a finite limit: 



Psui-v{Snj < 0, L) = W j- y = -— 
y p^nL p 



(AlO) 



The size of transverse fluctuations of the adsorbing 
walk is given by: 



kA5w.L)= PL{x)xdx/P,^„{5^,L) , (All) 
Jo 



where 



PL(x)xdx 



pV^ 



1 + ^(217^ + l)exp(l7^)[l- $(- 



(A12) 



r 



In the symmetric limit (5„, — 0: 



lu-{Sn^ - 0,L) = 78^-^1/2 



(A13) 



For Sw > the transverse fluctuations in leading order 
are independent of L: 



ltr{Sw >0,L) = 



2a' 



(A14) 



of positive steps is followed by a drift process towards 
the average behavior. If the number of steps in the drift 
process is aL, where < a < 1, then 

I'la^'^S^ > 0, L) = aL6^ . (A16) 



while for < 0, when there is a drift of the walk from 
the wall the transverse fluctuations grow linearly with L: 



ltr{Sw < 0, L) = SwL 



(A15) 



The maximal value of the transverse fluctuations 
^™"^(5u), L) for 5 < are in th e sam e ord er of magnitude 
as their average values in eqs(A13) and (A15). However 
for 6 > the maximal value is generally larger, then the 
average one in eq( |A14 ). In this case /[^"^(J^ > 0,L) is 
determined by a rare event, in which a large fluctuation 
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